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Abstract 

In order to maximize their fitness, individuals aim at choosing territories offering the most appropriate combination of 
resources. As population size fluctuates in time, the frequency of breeding territory occupancy reflects territory quality. We 
investigated the relationships between the frequency of territory occupancy (2002-2009) vs. habitat characteristics, prey 
abundance, reproductive success and parental traits in hoopoes Upupa epops L, with the objective to define proxies for the 
delineation of conservation priority areas. We predicted that the distribution of phenotypes is despotic and sought for 
phenotypic characteristics expressing dominance. Our findings support the hypothesis of a despotic distribution. Territory 
selection was non-random: frequently occupied territories were settled earlier in the season and yielded higher annual 
reproductive success, but the frequency of territory occupancy could not be related to any habitat characteristics. Males 
found in frequently occupied territories showed traits expressing dominance (i.e. larger body size and mass, and older age). 
In contrast, morphological traits of females were not related to the frequency of territory occupancy, suggesting that 
territory selection and maintenance were essentially a male's tasl<. Settlement time in spring, reproductive success achieved 
in a given territory, as well as phenotypic traits and age of male territory holders reflected territory quality, providing good 
proxies for assessing priority areas for conservation management. 
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Introduction 

Territory choice has a crucial impact on individual fitness since 
the acquisition of a highly suitable territory is essential for survival 
and successful reproduction [1,2]. Individuals should settle in the 
highest quality territories to maximize their fitness. High quality 
territories are those that contain an optimal combination of 
essential resources needed for reproduction and survival, such as 
nest-sites, food, and concealment from predators [3] . If individuals 
can move between territories, the ideal free distribution model 
(IFD) predicts that they should occupy habitat patches in 
proportion to the amount of available resources [4,5]. Thus, in 
areas of high habitat quality, territories are expected to be smaller 
than those in areas of lower quality habitats. Due to a roughly 
similar amount of resources available per territory this should 
result in equivalent reproductive success among territories across 
the area, insofar as an IFD model applies. 

The IFD model assumes that all individuals have the same 
competitive abilities, an assumption which is often violated [6]. 
The ideal despotic distribution model (IDD), on the other hand, 
takes the variations in competitive abilities among individuals into 
account and thus predicts that stronger, rather than weaker, 
individuals settle in higher quality territories [4,5] . Under an IDD 
scenario, reproductive success varies spatially, with greater success 
in territories of high quality. Both models postulate that the 
distribution of breeders in a heterogeneous environment is non- 
random in space and time: high quality territories tend to be 



occupied when population density is low, whereas low quality 
territories are only occupied during high population density. Thus, 
the number of times a territory is occupied over a given period of 
time (i.e. its occupancy rate) can be used as a reliable measure for 
territory quality [1,7-9]. 

Using territory occupancy frequency as a measure of territory 
quality provides a means of identifying the key habitat factors 
determining quality [1,10]. These factors can include any 
resources such as food supply, nesting sites and/or structural 
habitat variables. The identification of these key factors can 
constitute an essential step in species' conservation and manage- 
ment. 

Migratory bird species have to select a breeding territory every 
year. Typically, older individuals arrive first at the breeding 
grounds [11,12] and can therefore freely select their territory. In 
addition, early arriving individuals are often in better physiological 
condition than those arriving later [13]. The outcome of this 
settlement process is that the best territories are typically occupied 
by the most dominant, i.e. highest quality, individuals [14]. 

The competitive ability of an individual is difficult to assess, and 
usually phenotypic traits, such as age, body size or mass are 
considered. Older individuals often dominate younger individuals, 
since they are more experienced [2]. The same holds for larger 
and/or heavier individuals that are physically stronger [3,15]. 
Since individuals with high competitive abilities (dominant 
individuals) choose territories of higher quality, the frequency of 
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territory occupancy often correlates positively with phenotypic 
traits expressing the dominance of territory occupants [14]. 

Few studies have investigated the links between territory 
occupancy, breeding success, environmental characteristics, and 
parental phenotypic traits [16]. We studied the correlates and 
determinants of territory selection of hoopoes (Upupa epops L.) in 
southwestern Switzerland. Using the frequency of territory 
occupancy during 8 years as an indication of territory quality, 
we studied the relationships between territory occupancy, habitat 
parameters within territories, and individual characteristics of 
territory occupants at two spatial scales, and then evaluated the 
consequences of territory occupancy and individual qualities on 
reproductive success. Previous studies of the same population have 
shown that hoopoes mainly feed on mole crickets [Gryllotalpa 
gryllotalpa L.) [17,18], which are sought for amongst sparse ground 
vegetation [19]. We therefore predicted a positive relationship 
between frequency of territory occupancy vs. mole cricket 
occurrence and amount of bare soil, respectively. Secondly, we 
tested whether phenotypic traits of territory occupants were 
related to frequency of territory occupancy, predicting that larger, 
heavier and older individuals occupy territories that typically have 
a higher occupancy rate. Finally, we tested whether reproductive 
success was related to territory occupancy or to the phenotypic 
traits of territory occupants, in both cases anticipating positive 
relationships. Our ultimate goal was to establish more straightfor- 
ward cues to assess territory quality other than occupancy and thus 
to spatially identify priority areas for efficient conservation 
management. 

Methods 

Study species 

Hoopoes preferably inhabit semi-open, dry and sunny areas of 
southern Europe, north-western Africa and central Asia. Typical 
breeding habitats often include traditionally cultivated areas in 
central Europe [18], but high-intensity farmland can also be 
inhabited provided that essential resources are available [17,19]. 
Hoopoes are eclectic secondary cavity breeders, occupying hollow 
trees and walls, as well as nest boxes. Hoopoes often raise two 
broods a year [17]. After World War II hoopoes have undergone 
large declines in central Europe, including Switzerland. The 
species is considcr(-d to lie \uln(-rable in Switzerland [20], 
although populations have been increasing recendy [21]. 

Study area 

The study was conducted on the plain of the upper Rhone 
valley (Valais, 46.2°N, 7.4°E; 480 m above sea level) from April to 
August 2002-2009. The study area covered 62 km^, dominated by 
industrial farming, mainly consisting of fruit-tree plantations, 
vineyards and vegetable crops. For a more detailed description of 
the region see [21]. In the study area, 690 nest boxes designed for 
hoopoes (dimension 20x20x30 cm; entrance hole diameter: 
55 mm) were installed gradually from 1997 to winter 2002 at 
367 locations, with usually two nest boxes installed in a given 
building. The focal hoopoe population responded rapidly to the 
installation of nest boxes, with a ca. 2-fold increase in population 
size within a few years during the study period [2 1] . 

Sampling design 

Ethics statement. AU sites were visited with landowner 
permission. AU research protocols involving experimental activities 

with birds were approved by the Swiss Federal Office for the 
Environment (FOEN) and the Swiss Ornithological Institute and 
comply with the Swiss legislation. 



Territory occupancy. To determine territory occupancy we 
inspected nest boxes every second week throughout the breeding 
season. Ever)' year the inspections started before the first nest 
boxes were expected to be occupied (i.e. end of March). All 
occupied nest boxes were thereafter checked every third day to 
collect information about clutch size, number of fledglings, and 
phenology. For each of the 300 locations we added up the number 
of years in which they were occupied by breeding hoopoes, 
hereafter designated 'territory occupancy' (range 0-8 years, from 
2002-2009). We defined a territory to be occupied in a year if at 
least mic egg was laid in it in that year. Doubk" occupancies of a 
territory within a given year were not considered because, usually, 
at the time when second broods are initiated, many nest boxes are 
already occupied and consequently territory choice is already 
much restricted. 

Individual phenotypic traits. The majority of breeding 
hoopoes were caught using mist-nets or spring traps placed at the 
nest box entrance; females were also taken by hand from the nest 
box while brooding. Capture attempts started 4 days after the 
hatching of the first chick; they ended after fledging or after a 
maximum effort of 12.5 h for capture. The same standardized 
design was used every year. We determined sex and age (two 
classes: one year old and older, recognizable from moult patterns 
of wing feathers; Mosimann-Kampe unpublished data), and 
measured the lengths of bill, crest, the fifth primary feather (P5), 
the first tail feather (Rl) and tarsus, as well as body mass. AU 
captured individuals as weU as aJl nesdings were ringed [22]. 

Habitat mapping. From a total of 367 available breeding 
locations, only those that were occupied at least once by a 
breeding hoopoe between 2002 and 2008 were retained (n= 172). 
100 of them were chosen randomly, the habitat of which was 
mapped in 2009. A circle with a 300 m radius (corresponding 
approximately to a mean home range size of 40 ha [23]) was 
drawn around each of the 100 nest box locations. This circular 
area was considered to be a territory. Various habitat variables (see 
Table 1 and Table S 1 for descriptive statistics) were recorded at 30 
locations randomly selected within this circular territory in order 
to describe habitat features. In addition, we used a soil 
penetrometer to measure soil density at aU samphng points. To 
correct for the seasonal changes in soil density, we measured the 
soil resistance once a week (n = 15) between April and August at a 
single location in the middle of the study area (46°13'03"N 
07°20'41"E). The chosen location showed average soil conditions 
comparable to the main conditions prevailing at most of the other 
sampling points. Five measures each were taken at every sampling 
occasion. In addition, we used ground water table maps retrieved 
from pubUcly available data (Departement des Transports de 
I'Equipement et de I'Environnement, Etat du Valais & Centre de 
Recherche sur I'Environnement Alpin - CREALP) to look at the 
link between mole cricket occurrence and ground water depth. 

Mole cricket occurrence. To assess the relationship be- 
tween mole cricket oc[:urrence and habitat factors we sampled 
detection/ non-detection data of mole cricket traces at 97 plots, 
along with abiotic and biotic habitat factors (depth of ground 
water table, soil density, soil type, vegetation cover) in 2009. These 
plots were randomly selected using ArcGIS 9.x [24] within the 
same area where all other habitat variables were collected (see 
above), but exclusively within fruit tree plantations, because this is 
the favoured foraging habitat of hoopoes [23]. In order to ensure 
sufficient contrast, we used a stratified sampling design: plots were 
chosen at random in areas known to have a high and low ground 
water table, respectively (10 in each). Another 10 plots were 
randomly selected within areas with very graveUy soil. The 
remaining 67 plots were chosen randomly from the remaining 
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Table 1. Description of habitat variables recorded at nnapping locations for modeling hoopoe territory occupancy and mole 
cricket occupancy. 





Parameter 


Levels 


Description 


Habitat type 


apple tree plantation 




apricot tree plantation 


pear tree plantation 


arable field 


vineyard 


grassland 


river bank 


wood 


non-tarred road 




unsuitable area 


building, tarred road, open water 


Vegetation cover 




continuous (to nearest 10%) 


Mowing 


yes or no 


regular mowing of the driving track {only for fruit tree plantations) 


Ground 


mowing 


management of the vegetation strip underneath plantation trees (only for fruit tree 


management 




plantations) 




herbicide application 


strip underneath plantation trees (only for fruit tree plantations) 




mechanical veg. removal 


trees (only for fruit tree plantations) 




no treatment 


plantations) 


Soil type 


1} silty soil with no-till-limited presence of sand 


characterisation of top soil layer 


2} silty soil with obvious presence of sand 




3} silty soil embedded in a matrix dominated by gravel, stones or 




pebbles 






4) sandy soil where large structures such as gravel and pebbles 




are absent 






5) sand embedded in a matrix dominated by gravel, stones or 






pebbles 




6} all kind of humus rich soil (decomposed litter) 


Soil density 




continuous (0-15 in steps of 0.5) Five measures each at every sampling occasion 






using a soil penetrometer 



doi:l 0.1 371 /journal.pone.0097679.t001 



areas that showed average soil type and ground water table 
conditions. The presence of mole crickets was determined by 
searching, during 10 min, for underground galleries and holes at 
the soil surface along three 1 0 m transects separated by three tree 
rows. If at least one gallery or entrance hole was found, mole 
cricket presence was considered as detected at that visit. All plots 
were sampled four times in June, during the peak period of mole 
cricket activity [25]. We recorded habitat covariates we thought 
were relevant for mole cricket occurrence (Table 1), but also took 
into account covariates that might have affected detectability, 
namely vegetation height and weather variables (precipitation on 
the day preceding a visit; daily average temperature; daily 
maximum temperature, and daily duration of sunshine; all 
obtained from the meteorological station at Sion: 46°13'06"N, 
07°19'48"E). 

Statistical analyses 

Random territory choice and effects on productivity. All 

statistical analyses, except estimations of mole cricket occupancy 
(see below), were conducted using the software R [26]. To test 
whether hoopoes selected breeding locations randomly, we 
compared the frequency distribution of the observed occupancy 
of all territories that were occupied at least once from 2002-2009 
(n= 192) with the frequency distribution of occupancy that would 



be predicted in a scenario of random selection of territories. We 
used a simulation approach to determine the frequency distribu- 
tion of occupancy under random territory selection that takes into 
account the variable number of breeding pairs. Specifically, we 
randomly selected for each year from the list of the available 
territories the number of territories that corresponds to tiie 
number of breeding pairs in a given year and assigned them as 
occupied. We then tabulated occupancy of each territory, derived 
the frequency distribution and repeated this procedure 1000 times. 
The mean frequency distribution was then taken as the reference 
for random territory selection. We compared the two frequency 
distributions with the Pearson's chi-square test. 

To study the relationships between breeding phenology and 
reproductive success vs. territory occupancy we used linear mixed 
effects models. We used the hatching date of chicks as a proxy for 
the date of territory settlement and then modelled it using territory 
identity as a random effect, and the year of territory occupancy as 
a frxed effect. The inclusion of territory identity enabled us to 
include modelling of the hatching dates for every year, thus 
avoiding pseudo-replication. Reproductive success represented by 
the total number of fledglings per territory and year was modelled 
as a function of territory occupancy (i.e. number of times the 
territory was occupied between 2002 and 2009) using territory 
identity as a random effect and assuming a Poisson error 



PLOS ONE I www.plosone.org 



3 



May 2014 | Volume 9 | Issue 5 | e97679 



Breeding Territory Occupancy in Hoopoes 



100 



80 - 



CO 
<D 



•5 



60 



0) 

I 40 



20 




■ expected 
□ observed 



Number of years occupied 

Figure 1. Observed and expected frequency of territory occupancy. The comparison of observed and expected (according to a random 
selection scenario) frequency of territory occupancy of hoopoes (2002-2009) shows the deviation from a random territory selection pattern 
(X^ = 117.85, df=7, P<0.001). 
doi:1 0.1 371 /journal.pone.0097679.g001 



distribution. Furthermore, we tested whether the probability of 
second broods was related to the occupancy of territories. For this 
purpose we modelled the binary variable indicating the presence 
or absence of a second brood in every year and territory with 
territory occupancy, using territory identity as a random effect. All 
these models were compared to corresponding models without the 
variable 'territory occupancy', and model averaged parameter 
estimates were calculated based on AIC weights (Akaike's 
Information Criterion (AIC) [27]). 

Individual phenotypic traits. To evaluate whether pheno- 
typic traits of territory occupants were related to territory 
occupancy, we fitted linear mixed effects models (normal error 
distribution) with different morphological traits (biU length, crest 
length, length of P5 and first tail feather Rl, tarsus length, body 
mass) as response variables, and territory occupancy as an 
explanatory variable. Again, territory identity was included as a 
random effect. These models were evaluated for males and females 
separately given the existence of a slight sexual cUmorphism (males: 
n = 626; females: n= 758). As before, we compared the models 
with corresponding models without the variable 'territory occu- 
pancy' and used model averaged parameter estimates based on 
AIC weights for inference. 

Secondly, we tested whether the age of territory occupants (two 
categories: one year old and older) was related to territory 
occupancy. Linear mixed effect models with a binomial error 
distribution and territory identity as a random effect were used to 
model the binary response variable age. Modelling was performed 
separately for males (n = 553) and females (n = 654). The resulting 
models were again compared with a model without the variable 
'territory occupancy', and model averaged parameter estimates 
based on AIC weights were used for inference. 

To test whether the number of fledglings was linked to the 
morphology and age of breeding individuals we used linear mixed 



effects models with a Poisson error distribution. The number of 
fledglings of a given breeding adult in a year was modelled using its 
different morphological traits and age as explanatory variables, 
and the individual identity as a random effect. Confounding 
factors such as the year and date of the current brood, and 
whether it was a first or second brood were included in all models. 
The various models fitted were predefined in a sensible way and 
included either morphological traits and age separately, single 
morphological traits in combination with age, all morphological 
traits plus age, or none of these effects. Modelling was again 
performed for males (n = 397) and females (n = 434) separately. 
We computed model-averaged parameter estimates based on AIC 
weights. 

Mole cricket occurrence. Mole cricket detection/non- 
detection data was analysed using occupancy models [28] with 
software MARK [29]. We used average daily temperature, 
maximum daily temperature, amount of rainfall, duration of 
sunshine, sampling event, vegetation cover, and vegetation height 
to model detection probability (p). To model occupancy proba- 
bility (*?) we used depth of the ground water table, soil density 
(corrected for date, see above), soil type and vegetation cover as 
well as the quadratic terms of ground water table depth, soil 
density and vegetation cover, and two interactions, namely soil 
type*ground water table and soU type*soil density. We performed 
modelling in four steps. In a first step we used the full model for 
occupancy including all quadratic terms and interactions and 
explored the effect of different combinations of covariates for 
detection probability. In steps two and three we investigated the 
importance of different quadratic terms and interactions, respec- 
tively, on occupancy. Finally we evaluated the remaining 
combinations of variables that were neither involved in a quadratic 
term nor in an interaction on occupancy. The models were ranked 
in every step using AICc values (AIC corrected for small sample 
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Figure 2. Relationship between territory occupancy and 
individual characteristics. Relationship between territory occupancy 
(2002-2009) vs. P5 (fifth primary feather) length, tarsus length and body 
mass of territory occupants. IVIales: grey bars; females: white bars. The 
regression lines refer to males for which trends were statistically 
significant. 

doi:10.1371/journal.pone.0097679.g002 

size [27]), and the best models (AAICc ^ 2) were then selected for 
the next step. To get the best parameter estimates for the 
extrapolation of mole cricket occurrence across the study area we 
performed model averaging, subsequently using the best models of 
the final step, i.e. those models for which the sum of AICc weights 
was > 0.9. 

Relation of territory occupancy to habitat 
characteristics. To correct for seasonal changes of soil density 
we fitted a linear mixed effects model with soil density of the 
reference location (see above) as a response variable, date as an 
explanatory fixed variable, and visit as a random effect. This 
model revealed that soil density increased during the course of the 
season. Using the parameter estimates of this model we corrected 
all soil density measures for date. 

Using the recorded habitat variables, and extrapolated infor- 
mation about mole cricket occurrence (see above) we modelled 
territory occupancy (binomial data) of hoopoes in the randomly 
selected 100 territories using mixed effects models (territory 
identity number as a random effect, habitat variables, and the 
predicted probability of occurrence of mole crickets, as frxed). The 
inclusion of territory identity as a random effect accounted for 
repeated measurements at territories, as each territory was 
delineated by up to 30 sampUng points. We used the model- 
averaged parameter estimates obtained from the mole cricket 
occupancy modelling (see below), and point-specific habitat 
covariates to calculate the probability of mole cricket occurrence. 

We modelled territory occupancy in two main steps each using 
territory variables from two spatial scales (large scale = 300 m and 
small scale: a subset of the large scale = 200 m radius). First, we 
evaluated the effect of habitat type on territory occupancy. Since a 
habitat type could be assigned to all points, the sample size was 
3000 for the large spatial scale and 1388 for the smaller spatial 
scale. Secondly, we modelled structural habitat variables. Since 
they could not be assessed at all points, sample size was lower 
(large scale: n = 2378; small scale: n= 1153). Structural variables 
included soil density (corrected for date, see above), depth of 
ground water table, vegetation cover, and soil type as well as the 
probability of mole cricket occurrence as a measure of food 
abundance. We fitted a null model with an intercept only, a fuU 
model that included all structural variables and models that 
included only one structural variable. We also evaluated different 
combinations of structural variables and models with quadratic 
effects for soil density and vegetation cover. As fruit tree 
plantations are the main local foraging habitat of hoopoes [23], 
structural variables in that habitat may be the main determinants 
of overall habitat quality. Therefore, we tested the same 
combinations of models but this time only with structural variables 
measured in fruit tree plantations (large scale: n=1182; small 
scale: n = 634). 

Results 

Random territory choice and effects on productivity 

Hoopoes did not select territories at random (x'^= 117.85, 
df = 7, P<0.001, Fig. 1). If territory selection had been random, we 
would have expected a greater number of territories to be 
occupied one to three times and a lower number of territories to be 
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occupied more than four times than observed. The hatching date 
decreased noticeably with the frequency of territory occupancy 
(P=-2.44, SE = 0.39, AIC weight >0.99); thus, territories that 
were occupied more often were also occupied earlier in the season. 
The number of fledglings was also positively related with the 
frequency of territory occupancy (P = 0.07, SE = 0.01, territory 
random effect (variance) = 0.05, AIC weight >0.99); thus more 
frequently occupied territories were overall more productive than 
less frequently occupied territories. The probability of rearing a 
second brood increased with the frequency of territory occupancy 
(P = 0.40, SE = 0.08, territory random effect (variance) =0.97, 
AIC weight >0.99), thus, the larger number of fledglings that were 
raised annually in high quality territories was mostly the result of a 
larger annual number of successfully reared broods. 

Individual phenotypic traits 

The average age of breeders was strongly related to territory 
occupancy in males (P = 0.23, SE = 0.04, territory random effect 
(variance) <0.01, AIC weight >0.99), but not in females (P = 0.04, 
SE = 0.04, territory random effect (variance) <0.01, AIC weight 
= 0.65). The relationship was positive in both sexes, thus older 
individuals tended to settle in territories that were more frequendy 
occupied. 

There were positive relationships between territory occupancy 
and morphological characteristics in males, but not in females. In 
males, morphological traits increased with increasing territory 
occupancy, with three of them statistically supported: P5 length, 
tarsus length and body mass (Table 2). Thus, males present in 
territories of higher occupancy were larger and heavier than males 
present in rarely occupied territories (Fig. 2). By contrast, the null 
model was often supported in females, with only one morpholog- 
ical trait whose 95% confidence interval was non-overlapping with 
0 (Table 2). This trait (tarsus length), however, had a negative 
estimate, indicating that small sized females would have a better 
chance to occupy a territory. There was no link between individual 
reproductive success and morphological traits and/or age. None of 
the models using morphological traits and/or age to model 
fledgling numbers was statistically supported (Table S2). 

Mole cricket occurrence 

Mole cricket occupancy models with interactions (soil type*- 
ground water table and soil type*soil density) were weakly 
supported by the data (AAICc = 13.72 compared to overall best; 
Table 3). The quadratic terms of soil density and ground water 
table seemed to be important and were included in the best 15 
models. The averaged parameter estimates showed that the 
probability of detecting mole crickets decreased with increasing 
ground vegetation cover (P = —0.63; SE = 0.16), with a decreased 
duration of sunshine (P = 0.43; SE = 0.14), and varied slighdy 
among sampling events. 

The parameters most relevant for mole cricket occurrence were 
soil type, soil density, depth of ground water table and vegetation 
cover, as well as the same two above-mentioned quadratic terms 
(Fig. 3). It is striking that mole cricket occupancy was similar in all 
soil types except soil type 3 (sUty soil embedded in a matrix 
dominated by gravel, stones or pebbles), in which mole crickets 
had a considerably lower probability of occurrence. Soil type 6 
(soil rich in humus) did not occur in any sampled fruit tree 
plantation and was therefore not included. Mole cricket occur- 
rence declined with increasing depth of the ground water table 
(Fig. 3). The relationship with soil density was more complex; mole 
cricket occurrence was higher in soft and hard soils compared to 
medium soils (Fig. 3). This was most pronounced with soil type 3. 
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Table 3. Mole cricket occupancy models. 



Model 


AAlCc 


Wi 


K 


Deviance 


f (soilt+dens+gw+dens^) p (veg+sun) 


0.00 


0.46 


11 


401.15 


^ (dens+gw+dens^) p (t+veg+sun) 


3.35 


0.09 


10 


407.05 


4^ (dens+dens^) p (t+veg+sun) 


4.51 


0.05 


9 


410.72 


f (soilt+dens+dens^) p (t+veg+sun) 


4.79 


0.04 


13 


400.62 


*P (soilt+dens+veg+dens^) p (t+veg+sun) 


5.01 


0.04 


14 


398.09 


*P (dens-veg+gw+dens^) p (t+veg+sun) 


5.06 


0.04 


11 


406.21 


^ (dens+gw+dens^+gw^) p (t+veg+sun) 


5.11 


0.04 


11 


406.26 


*P (dens+gw+dens^) p (veg+sun) 


5.51 


0.03 


7 


416.55 


4^ (dens+veg+dens^) p (t+veg+sun) 


5.86 


0.02 


10 


409.57 


*P (soilt+dens+dens^) p (veg+sun) 


6.31 


0.02 


10 


410.02 


4^ (soilt+dens+veg+dens^) p (veg+sun) 


6.45 


0.02 


11 


407.60 


4^ (soilt+dens+gw+dens^+gw^) p (t+veg+sun) 


6.46 


0.02 


15 


396.71 


4* (soilt+dens+veg+gw+dens^+gw^) p (t+veg+sun) 


6.66 


0.02 


16 


394.01 


4* (soilt+dens+gw+dens^) p (t+veg+sun) 


6.72 


0.02 


14 


399.80 


4* (dens+dens^) p (veg+sun) 


6.75 


0.02 


6 


420.12 


4^ (dens+veg+gw+dens^+gw^) p (t+veg+sun) 


6.82 


0.02 


12 


405.34 


4* (dens+veg+gw+dens^) p (veg+sun) 


7.07 


0.01 


8 


415.73 


4^ (dens+gw+dens^+gw^) p (veg+sun) 


7.08 


0.01 


8 


415.74 


4^ (soilt+dens+veg+gw+dens^) p (t+veg+sun) 


7.22 


0.01 


15 


397.47 


4^ (soilt+dens+gw+dens^+gw^) p (veg+sun) 


7.56 


0.01 


12 


406.08 


4^ (soilt+dens+veg+gw+dens^+gw^) p (veg+sun) 


7.67 


0.01 


13 


403.51 


4^ (dens+veg+dens^) p (veg+sun) 


8.00 


0.01 


7 


419.04 


4^ (soilt+dens+veg+gw+dens^) p (veg+sun) 


8.47 


0.01 


12 


406.99 


4* (dens+veg+gw+dens^+gw^) p (veg+sun) 


8.64 


0.01 


9 


414.85 



Model selection summary for mole cricket occupancy (4*) and detection probability (p) in response to habitat parameters. Shown are the differences between the best 
and the current model (AAlCc), the AlC weight of the current model (iVj), the number of estimated parameters (K) and the deviance. The best 1 5 models {J^Wi = 0.9) are 
shown on top separated from the others by a horizontal line. 

Covariates: dens = soil density, gw = height of ground water table, soilt = soil type, sun = daily sunshine duration, t= sampling occasion, veg = vegetation cover, 
dens^ = quadratic term for soil density, gw^ = quadratic term for ground water table. 
doi:l 0.1 371/journal.pone.0097679.t003 



As expected, vegetation cover had only a weak impact on mole 
cricket occupancy (Fig. 3). 

Relation of territory occupancy and habitat 
characteristics 

At both spatial scales, the simplest model including none of the 
habitat variables was the best, suggesting that the recorded habitat 
characteristics and mole cricket occurrence were not related to 
territory occupancy (Table 4). This was also true when only fruit 
tree plantations were considered (Table S3). 

Discussion 

This study shows that hoopoes select their territories in a non- 
random manner, indicating spatial heterogeneity regarding the 
availability of crucial resources. In agreement with our expecta- 
tions, territory settlement occurred earlier in the season and 
reproductive success was higher in frequently occupied territories. 
However, our models were unable to identily the environmental 
determinants of territory occupancy: neither habitat type, nor 
vegetation characteristics, nor mole cricket occurrence correlated 
with the frequency of territory occupancy. This may be due, at 



least in part, to some artifact linked to the method used; we shall 
return to this point later on. 

In line with the predictions of the ideal despotic distribution 
model (IDD), male hoopoes breeding in often occupied territories 
were larger and heavier than birds breeding in territories that were 
rarely occupied. However, no such relationship was apparent in 
females. The demonstration that frequently occupied territories 
had a higher reproductive output corroborates the findings from 
other studies [1,3,7-9,30-32]. Frequendy occupied territories were 
more likely to support second broods, which was key to increased 
territory-specific reproductive success. The difference in repro- 
ductive output among territories can either originate from 
differences in territory quality per se, or from differences in the 
quality of territory occupants [33]. Since hoopoes are short-lived 
(annual survival probability is about 40%, [22]), difiFerent 
individuals have contributed to the high breeding success in the 
territories of high occupancy. This suggests that differences in 
territory quality were of importance. However, experiments are 
necessary to disentangle effects of territory and individual quality. 
The non-random selection of breeding territories with a clear role 
of male phenotypic traits suggests that a hierarchical system is 
operating. This pattern again conforms to the IDD model. 
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Figure 3. Occupancy probability of mole crickets in response to 
key habitat properties. Occupancy probability of mole crickets in 
response to a) depth of the ground water table, b) soil density, and c) 
surface vegetation cover. Model-averaged predictions are shown for 
five different soil types occurring at the sampling sites (Table 1). Closed 
circles represent soil type 1, open circles represent soil type 2, open 
squares represent soil type 3, open diamonds represent soil type 4 and 
open triangles represent soil type 5 (see Table 1 for the description of 
soil types). 

doi:10.1371/journal.pone.0097679.g003 



The relationships between territory occupancy and individual 
characteristics (age, size, body mass) difTered between the genders. 
This suggests that males and females may rely on different criteria 
for operating territory selection and deciding about reproductive 
strategy. In multiple-brooding bird species it is commonly the task 
of the male to judge territory quality and try to monopolize the 
best available sites [1 1], while females evaluate the quality of males 
rather than the quality of territories [30]. The patterns observed in 
hoopoes seem to match this view. First, it is exclusively the task of 
males to secure favorable territories, as commonly observed in 
other species [11]. Males are also more phUopatric to their 
breeding territory than females [34] . Altogether, this suggests that 
dominant males monopolize the best territories and avoid 
dispersal. 

There are two mechanisms that could explain why frequendy 
occupied territories were occupied earlier in the season by older 
and larger hoopoe males: first, older birds are usually more 
experienced and thus have a higher hierarchical status enabling 
them to outcompete younger individuals [2,3,7,12,15,31,35]. 
Secondly, older migrating birds are more experienced, usually 
arriving earlier on the breeding grounds than younger ones, and 
are therefore free to choose the best territories [2,11,12]. 

Hoopoe males in frequently occupied territories were larger and 
heavier than males present in rarely occupied territories. Theory 
predicts that dominant individuals will occupy high quality 
territories [4] , and thus we can conclude that body size and body 
mass are important for defining the hierarchical status of hoopoe 
males. It is, however, not clear whether body size is direcdy 
beneficial in terms of agonistic confhcts for territory acquisition, or 
indirectly as a result of the earlier arrival date [2,15]. First-arrived 
males may be more efficient migrators that benefit from 
precedence in territory settlement. In contrast, body size traits of 
females were not positively correlated with territory occupancy. As 
females lose body mass during incubation [36], a possibly existing 
relationship may have been blurred, since not all females could be 
captured at the same nesting stage. For the other female traits, 
however, analogous spurious findings are unlikely. 

Contrary to our expectations, we did not find any connection 
between frequency of territory occupancy and habitat character- 
istics or the occurrence of the main prey, mole crickets. This held 
true both when considering two spatial scales, and when only the 
favorite foraging habitat type (fruit tree plantations) was consid- 
ered. These results could be interpreted as if all territories were 
equal regarding habitat and prey availability. Indeed in monot- 
onous human-made ecosystems like fruit tree plantations, the 
variation of some habitat features is small due to similar 
management practices (Table SI) and the quality of the nesting 
site per se might be more important than the seemingly 
homogenous habitat around the nesting sites. However, as 
foraging hoopoes showed a clear preference for certain habitat 
types [23], we rather suggest that it is more likely that extant 
differences were simply not detected via our approach. First, it can 
be argued that we did not consider the most relevant habitat 
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variables. This is quite unlikely, yet, given that our habitat 
descriptors were chosen from fine-grained, radio-tracking infor- 
mation about habitat selection [23]. Neither can changes in the 
farmland matrix during the course of the study be inferred, 
because fruit tree plantations, predominant in the study area, have 
a low renewing rate. Secondly, a recent study has showed that the 
proportion of mole crickets in chick's diet actually increases with 
territory occupancy (Guillod, unpublished data), suggesting that 
the abundance of the locally most profitable prey does link to 
territory occupancy and is thus a prime determinant of territory 
quality. We might simply have failed to evidence such a link 
because we could not measure mole cricket abundance directiy 
and had instead to use mole crit:ket occurrence as a measure of 
food availability, which probably does not well reflect local 
abundance. Furthermore, we predicted mole cricket occurrence 
indirectly, by using habitat variables, and any such extant link may 
have been blurred by too mui:h variation around the signal. 
Finally, and perhaps most importantly, the true distriljution 
pattern of mole crickets within hoopoe territories is likely to be 
very patchy, imposing additional difiiculty for a proper estimation 
of both spatial distribution and local abundance. In efiect, radio- 
tracking of foraging hoopoes indicate the existence of small areas 
where prey abundance is very high, forming mole cricket clusters 
[23]. Finally, our assumption that hoopoe territories are circular, 
with the breeding location at tlu; c(;nt(;r, is an extreme over- 
simpUfication. A better, but logisticaUy challenging approach 
would be to determine the true boundaries of any occupied 
territories by means of radio-tracking and then to measure mole 
cricket abundance and habitat in the (mtire territory. 

The occurrence of mole crickets dc'pends on soil structure. 
These insects preferred small-grained soils with lots of sand rather 
than gravelly soils and soils with limited amounts of sand. They 
also showed some preference for wet/humid soUs. Soil density also 
influences the occurrence of mole crickets, but the pattern is 
difficult to interpret. We expected mole crickets to prefer soils of 
medium density because a soft soil might cause mole cricket 
galleries to txjUapse, while in hard soils mole crickets cannot 
borrow galleries. Yet, we observed a high probability of 
occurrence in very soft and hard soils. Mole cricket occurrence 
further declined with increasing vegetation cover. This may be due 
to a preference of mole crickets for warm soils to ensure an optimal 
development of their clutches: short, sparse- vc-getation allows 
better warming of the soil than dense, high vegetation. 

Hoopoe territories clearly difiered in quality, as inferred from the 
non-random pattern of occupancy, relative reproductive output, 
and from the despotic distribution of male phenotypes. The pattern 
of territory occupancy and male hierarchical status estimated from 
phenotypic traits provide spatially-explicit information that might 
be important for setting conservation priorities: conservation action 
must focus on those areas in which territories are more frequently 
occupied (Fig. SI) and that are inhabited by dominant males. These 
frequentiy occupied territories produced most of the offspring, and 
as productivity is an important driver of hoopoe population 
dynamics [22], these territories are especially important for the 
survival of this population. Although occupancy data might be easily 
retrieved, they require a lengthy time span, which may represent a 
serious handicap in conservation projects. ReKance on more 
straightforward surrogates of territory occupancy represents a good 
alternative. As demonstrated here for the hoopoe, useful immediate 
proxies would be the phenology of territory settlement, the 
reproductive success achieved in a given territory, as well as the 
hierarchical status of a male territory holder. 



Supporting Information 

Figure SI Spatial distribution of territory occupancy. 

Map of the study area showing the different territories (black dots). 
The diameter of the dots corresponds to the number of years a 
territory was occupied from 2002-2009. For illustration purposes 
also territories that were never occupied from 2002-2009 are 
shown (red dots). The Kernel Density Tool of ArcMap 10.1 was 
used to interpolate the occupancy pattern over the study area. 
Areas where territories of high occupancy are aggregated are 
highlighted in white, areas with low occupancy in dark grey. Four 
large-scale high quality areas can be distinguished. 

gpG) 

Table SI Descriptive statistics of the continuous habitat 

variables. Basic statistical parameters of the different continuous 
habitat variables considered as predictors for territory quality. 
Shown are arithmetic means (Mean), standard errors (SE), 
minimum, maximum and the coefficient of variation (CV) for 
the radii of 300 m and 200 m. As habitat variables were recorded 
at 30 sampling points per territory, the statistical parameters 
shown in the table were calculated from the territory means that 
were calculated in a first step. 
(DOCX) 

Table S2 Influence of parental individual characteris- 
tics on breeding success. Model selection summary of the 
effects of 1) different morphological traits and 2) different 
morphological traits plus age on individual breeding success (i.e. 
number of fledglings) of male and female hoopoes, evaluated by 
linear mixed effects models. Besides the morphological traits, all 
models include the following co-variates: year, egg-laying date (of 
the current brood) and whether the current brood was a first or a 
second brood. Shown are the differences between the best and the 
current model (AAIC), the AIC weight of the current model (w;), 
the number of estimated parameters {E), the model deviance and 
the territory random effect variance (a^). 
(DOCX) 

Table S3 Environmental determinants of territory 
occupancy within tree plantations only. Model selection 
summary of the effects of habitat variables, in fruit tree plantations 

only, on the frequency of territory occupancy within a radius of 
300 m (n = 1 182) and 200 m (n = 634). Shown are the differences 
between the best and the current model (AAIC), the AIC weight of 
the current model (;»;), the number of estimated parameters (A), 
the model deviance and the territory random effect variance (ct^). 
(DOCX) 
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